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When a thin elastic sheet is confined to a region much smaller than its size the morphology of the 
resulting crumpled membrane is a network of straight ridges or folds that meet at sharp vertices. 
A virial theorem predicts the ratio of the total bending and stretching energies of a ridge. Small 
strains and curvatures persist far away from the ridge. We discuss several kinds of perturbations 
that distinguish a ridge in a crumpled sheet from an isolated ridge studied earlier (A. E. Lobkovsky, 
Phys. Rev. E. 53 3750 (1996)). Linear response as well as buckling properties are investigated. We 
find that quite generally, the energy of a ridge can change by no more than a finite fraction before 
it buckles. 
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I. INTRODUCTION 

There is an abundance of phenomena involving strong 
deformation of thin elastic membranes that span a wide 
range of scales. On the microscopic scale, phospholipid 
membranes behave like a solid below a two-dimensional 
(2D) freezing point Q. Some inorganic compounds such 
as graphite oxide || and molybdenum disulphite || also 
behave like elastic membranes on scales large compared 
to the interatomic spacing. The graphite oxide sheets 
can be collapsed in solution by inducing an effective at- 
tractive interaction between distant parts of the sheet. 
Molybdenum disuphite has also been observed in a "rag" 
phase that looks similar to a crumpled piece of paper. 
Mechanical properties of macroscopic thin elastic plates 
and shells undergoing large deformations are important 
in engineering of safety structures and packaging ma- 
terial development 0. 

Stability and post-bucking properties of thin shells and 
plates have been a subject of intense investigation Q. 
Few general results have been derived however. Due to 
the complexity of the equations which describe large de- 
flections of thin plates (two quadratic fourth order par- 
tial differential equations), rigorous proofs are difficult 
to achieve. Difficulties in treating thin shells and plates 
arise due to the fact that a small parameter related to 
the thickness of the shell multiplies the highest derivative 
term in the equations JtJ . It is well-known that this fact 
gives rise to a variety of boundary layer phenomena in 
the bending of thin shells Many different types of 
boundary conditions that lead to a boundary layer have 
been analyzed. They include bending moments shear 
forces pp| ] and free boundary conditions with distributed 
bending moments [ pd| . A common feature that emerges 



from these studies is that membrane stresses become con- 
fined to the boundary layer region whose size vanishes as 
some power of the shell thickness. 

It has been suggested recently by the author and oth- 
ers [ fl2"|Jlg[| that membrane stresses in a strongly crumpled 
elastic sheet are also confined to a set of straight ridge sin- 
gularities or folds. These ridge singularities, which were 
shown to arise under quite general conditions, constitute 
another case of the boundary layer phenomena in thin 
plates. A scaling law for the ridge width as a function 
of the plate thickness had been established with the use 
of an energy scaling argument | ^] and a boundary layer 
analysis of von Karman thin plate equations O]. Elas- 
tic energy was found to be confined in the ridges and to 
scale as 1 /3 power of the size of the ridge for a fixed plate 
thickness. Other scaling laws such as the dependence of 
the ridge width and energy on the dihedral angle were 
also investigated. 

In this article we extend the analytical and numeri- 
cal study of the ridges to explore whether the results 
obtained for isolated ridges can be successfully applied 
to a network of interacting ridges in a crumpled elastic 
sheet. In Sec. II we introduce the concept of the "mini- 
mal ridge" that refers to the necessary and sufficient con- 
ditions leading to the formation of the ridge singularity 
in the limit of the vanishing plate thickness. Deviations 
from these "minimal" boundary conditions in a crumpled 
sheet can then be treated, at least in the first approxima- 
tion, as perturbations to the "minimal" ridge. In Sec. Ill 
we discuss a "virial theorem" for ridges that is a direct 
consequence of the energy scaling argument. It provides 
a useful tool in testing the confinement of elastic energy 
and the degree to which ridge interaction influences it. In 
Sec. IV we present a scaling argument and an extension 
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of the asymptotic analysis of the von Karman equations 
that establish that there is a long range decay of the 
longitudinal ridge curvature and the transverse stress far 
away from the ridge. The energy in these ridge "echoes" 
is negligible compared to the energy of the ridge. A per- 
turbation framework for dealing with external forcing of 
the minimal ridge is presented in Sec. V. A particular 
case of the ridge compression is studied in detail and the 
scaling of the elastic energy correction is found. We show 
in Sec. VI that the small echo strains and curvatures are 
sufficient to change the energy of a far away ridge by a 
finite fraction. Numeric evidence corroborating some of 
the claims made in the preceeding chapters is presented 
in Sec. VII. Finally, implications of the ridge properties 
for the crumpling problem as well as future work are dis- 
cussed in Sec. VIII. 



II. THE "MINIMAL" RIDGE 

To facilitate a study of the ridge singularity one must 
first write down equations that describe large deflections 
of thin plates. Second, an understanding of the condi- 
tions that lead to the formation of the ridge singularity 
must be achieved. It was suggested in Ref. 1 13 that the 



existence of the sharp vertices where both radii of cur- 
vature are of the order of the sheet thickness is a neces- 
sary and sufficient condition for the formation of ridges. 
Ridges connect these points of high curvature. A min- 
imal way to create a ridge, therefore, Ref. Jl3| argued, 
is to introduce sharp points at the boundary of a flat 
rectangular piece of elastic material by requiring that its 
boundary follow a frame that has a sharp bend. The 
von Karman equations that describe large deflections of 
thin clastic plates can be used to deduce the asymptotic 
behavior of the ridge solution in this simple geometry. 

Let us recall the boundary value problem that ex- 
hibits the ridge singularity Consider a strip made 
of isotropic homogeneous elastic material with Young's 
modulus Y and Poisson ratio v. It has thickness h and 
width X. The points are labeled by the material co- 
ordinates (x,y) so that x E (— X/2, X/2). The strip ex- 
tends in the y-direction. Normal forces arc applied to the 
long boundaries x — ±X/2 in such a way as to force the 
boundary to follow a rectilinear frame that has a sharp 
bend at y = 0. The bend dihedral angle is ir — 2a (so 
that a = corresponds to a flat strip) . The shape of the 
strip and the elastic stresses are found from a solution to 
the non-dimensionalized von Karman equations 



V 4 /=[x,/] 



A 2 V 4 X = --[/,/] 



(la) 



(lb) 



Here V 4 = V 2 V 2 and a square bracket [a, b] denotes a 
symmetric contraction of the second derivatives of the 
fields a and b 



[a,b] = e af _ l e l3l/ (d a d l 3a)(d ll d v b) 

_ d 2 a d 2 b d 2 a d 2 b ^ d 2 a d 2 b 
dx 2 dy 2 dy 2 dx 2 dxdy dxdy 



(2) 



Here e a p is the antisymmetric two by two tensor. Sum- 
mation over repeated indices is implied. 

All lengths are measured in terms of the strip width 
X and energies in terms of the bending rigidity k = 
Yh 3 1 (12(1 — v 2 )). The first von Karman equation is the 
statement of the local normal force equilibrium. The sec- 
ond one has a purely geometric origin. It simply states 
that Gaussian curvature — | [/, /] acts as a source for the 
stress field. The small parameter A ~ h/X is propor- 
tional to the dimensionless thickness of the sheet. Here 
f{x,y) and x(x,y) are the potentials whose derivatives 
give the curvatures C a p and the two dimensional in-plane 
stresses a a/ 3 via 



C a p = X- 2 d a df3f 



cr Q/3 = kX 2 tanZpudpduX- 



(3) 



(4) 



To clarify the meaning of the curvature tensor C a p we 
note that its eigenvalues are the inverses of the principal 
radii of curvature of the sheet. The sheets assumes a con- 
formation that minimizes the elastic energy consisting of 
the bending and the stretching parts (measured in the 
units of k) 



E s tr — A Z 



dxdy[V 2 ff 



dxdy[S7 2 X ? 



(5a) 



(5b) 



There are two ways to supply boundary conditions for 
the von Karman plate equations consistently. First, one 
could specify a Kirchoff type condition on the functions 
/ and x an d their derivatives in terms of the material 
coordinates. Second, a prescribed shape in the embed- 
ding space corresponds to clamping the boundaries of the 
strip or simply supporting them. Only the first type of 
the boundary conditions is tractable in general since the 
relation between the derivatives of the functions / and x 
and the shape of the sheet in the embedding space is non- 
linear. In addition, even if one succeeds in translating the 
boundary conditions that make reference to the embed- 
ding space into the language of / and x, they will change 
when the thickness of the sheet is varied or external forces 
are applied. Therefore, only when the stresses and cur- 
vatures are specified at the boundary can one make any 
analytical progress. We must remark at this point that 
the effective boundary conditions for a ridge in a crum- 
pled sheet are not of this type. Stresses in the facets are 
non-zero. More importantly, they depend on how each 
ridge is stressed by the rest of the sheet. 
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The "minimal ridge" boundary conditions are given in 
terms of the curvature at the boundary. The only re- 
quirement is that there be two sharp points on each long 
boundary of the strip. The stresses and torques vanish 
at the boundary x — ±X/2 (except at the singular point 
y = 0) 



d a dpx = 9 a dpf = 0. 



(6) 



The sharp vertices introduced at the boundary can be 
mathematically expressed as singularity in the curvature 
at the boundary points y = and x = ±X/2 



d 2 f 

= "<%)■ 



(7) 



The coefficient a in Eq. ([?]) is exactly the bending an- 
gle alpha of the frame which is equal to the half of the 
difference of the dihedral angle of the frame from n [|l3| . 

The motivation behind seeking the "minimal ridge," is 
that a number of the asymptotic properties of the ridge 
singularity arc independent of the details of the boundary 
conditions. These include the exponents in the asymp- 
totic thickness scaling of the ridge curvature, elastic en- 
ergy and other quantities. Dihedral angle scaling expo- 
nents are also independent of the details of the bound- 
ary conditions JDjfl . Other ridge properties, that do de- 
pend on the details of the boundary conditions (the lon- 
gitudinal stress supported by a ridge, for example), can 
then be found perturbatively, at least in the first approx- 
imation. The following two sections explore additional 
boundary conditions-independent asymptotic properties 
of the "minimal ridge" that are useful in the investiga- 
tion of the effects of the crumpled sheet environment on 
the "minimal ridge." 



G = Yh is the two dimensional stretching modulus of a 
sheet of thickness h made of elastic material with Young's 
modulus Y. For the purposes of an energy scaling ar- 
gument we may ignore the second terms in the expres- 
sions for the e nerg ies. The argument given by Witten 
and Li in Ref. JlJ] estimates these energies in terms of 
a characteristic curvature C, for example, the transverse 
curvature in the middle of the ridge. According to that 
argument, if the length of the ridge, i.e. the distance be- 
tween the vertices, is X, then the characteristic strain 
7 ~ (CX)~ 2 exists in the ridge region of width C _1 so 
that the total stretching energy in the ridge is approx- 
imately E stT ~ G-f 2 (X/C) ~ GX- 3 C- 5 . Similarly, the 
bending energy is given by the characteristic ridge cur- 
vature C, via E hcnd ~ kC 2 (X/C) ~ nXC. Ref. fjj then 
argued that C is the only parameter that characterizes 
the ridge. Also, since the bending and the stretching 
energies both vary as a power of C, they must be compa- 
rable when the total energy is minimized. It immediately 
followed that C ~ (k/G)- 1 / 3 ^- 2 / 3 - h^^X- 2 ' 3 . 

Another important conclusion of the energy scaling ar- 
gument emerges when we consider the derivative of the 
energy with respect to the parameter C . It must vanish 
for the value of C at which the minimum total energy is 
achieved. Since the energies depend on powers of C we 
obtain the following statement 



dE dE str , dEbend 



-5E, 



dC dC dC 



str -^bcnd 

~c 



= o, 



(9) 



III. VIRIAL THEOREM 

A study of the "minimal ridge" [jl3) revealed that as 
the thickness of the sheet vanished, the shape developed 
a sharp crease. Details of how that singular limit is ap- 
proached were found. In particular, the elastic energy 
concentrates in a small region around the ridge of the 
size given by the characteristic ridge curvature C. This 
fact, allows one to construct an energy scaling argument 
that yields the asymptotic scaling behavior of the ridge 
curvature C and the elastic energy E Q . 

For the purpose of the energy scaling argument given 
in Ref. |l4|] , let us rewrite the expressions for the bending 
and the stretching energies Eqs. (^|) in terms of the princi- 
pal strains 71 and 72 and the principal radii of curvature 
Ri = 1/Ci and R 2 = 1/C 2 . We get £jj] 

£bc„d = ~ J dS (/s(Ci + G 2 ) 2 + k g CiC 2 ) (8a) 
E stI = \ f dS (G( 7 i + l2 f + G s7l72 ) . (8b) 



which leads to a "virial theorem" for ridges -Ebend = 5^str 
analogous to the virial theorem relating kinetic to poten- 
tial energy in celestial mechanics fig . The validity of 
this virial theorem rests on rather general grounds. The 
only requirements are that a) the energy be rigorously 
expressible a sum of a bending and a stretching contri- 
butions and b) each of these varies as a power of a free 
parameter, here C. Thus a failure of the virial theorem 
would indicate that the energies did not vary with the an- 
ticipated powers of C. Such a failure would be expected 
if G -1 were comparable to other lengths in the problem, 
such as the size of the sheet X. 

We have performed a numerical test of the virial theo- 
rem using a lattice model of an elastic sheet after Seung 
and Nelson |l7]]. A numerical verification of the virial 
theorem for the tetrahedral shape is presented in Sec. 
VI, here we only remark that the agreement is better 
than few percent. If a view of the crumpled membrane 
as a network of ridges that meet at vertices is correct, 
the tetrahedron virial theorem suggests that ridges in a 
crumpled sheet are well defined objects to which a scaling 
argument can be applied. 
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IV. LARGE DISTANCE BEHAVIOR OF THE 
RIDGE SOLUTION 



To ascertain how much ridges influence one another in 
a crumpled sheet one must determine how the boundary 
layer solution, that is valid in progressively narrow region 
around the ridge midline, joins onto the flat solution far 
away from the ridge. The following development is moti- 
vated by an observation that in a numerical implementa- 
tion of a ridge, the "echo" strains and curvatures persist 
far away from the ridge, despite the fact that most of the 
energy is concentrated into a region of asymptotically 
vanishing width. In fact, the characteristic decay length 
L of the "echo" strains and curvatures may diverge in the 
A — > limit. In this section we present an energy scal- 
ing argument for the existence and asymptotic scaling of 
the decay length L for the "minimal ridge" configuration. 
This conclusion is put on a more rigorous footing by an 
extension of the asymptotic analysis of the von Karman 
equations that includes the matching condition between 
the boundary layer and the large distance solution. The 
scaling of the long distance solution is later supported by 
numerical evidence. 

We imagine cutting the ridge along its midline into two 
identical parts. For the boundary value problem defined 
on one of resulting semi-infinite strips, the ridge effec- 
tively introduces some complicated boundary conditions 
applied at the cut. The detailed form of these boundary 
conditions is not important for our purposes. It is signif- 
icant, however, that they have no singularities. For the 
purposes of an energy scaling argument, we imagine that 
the short side of the semi-infinite strip is slightly bent so 
that its middle is displaced in the normal direction by an 
amount a as shown on Fig. 1. 




FIG. 1. Gentle curvature at a short end of a semi- 
infinite strip decays at a characteristic distance L. 

This curvature decays to zero at a characteristic dis- 
tance L from the short end of the strip. This length is 
set by the competition of the bending energy that favors 
quick decay of the curvature and the stretching energy 
that is smallest when the decay length is large. The bend- 
ing energy is given by the integral of the squared mean 
curvature over the area of the strip. If the decay length 
L is much larger than the vertical displacement a, the 
dominant curvature is in the x (short) direction and is 
approximately C xx ~ a/X 2 . The bending energy is then 
-Ebond ~ kC xx XL ~ na 2 X~ 3 L. Here n is the bending 



rigidity and XL is the area of the strip where the de- 
formation exists. The strain created by the decay of the 
transverse curvature C xx is due to the fact that the mid- 
dle of the sheet is inclined by a small angle of order a/L 
and thus the projected length of the line that bisects the 
strip is shorter than the length of the same segment of 
the boundary. The length mismatch creates a character- 
istic strain in the y (long) direction ^ yy ~ a 2 /L 2 . Thus 
the stretching energy is E stl - ~ Gj 2 y XL ~ nh~ 2 a 4 XL~ 3 . 
Here we used the fact that the 2D stretching modulus G 
is related to the bending rigidity via G ~ n/ti 2 . Since 
both kinds of energies vary as a power of L, they must 
be comparable when the total energy Sbend + ^str is min- 
imized. We thus, obtain the scaling of the decay length 
L ~ X(afh) 1 / 2 . The displacement a is determined by the 
scaling properties of the ridge. The asymptotic scaling of 
C xx ~ a/X 2 must be consistent with that of the the lon- 
gitudinal ridge curvature d 2 f/dx 2 ~ A 1 / 3 jX. Therefore 
a ~ XX 1 / 3 which yields the scaling for the decay length 
L and the elastic energy in the ridge "wing" in the A — > 
limit 



L ~ XX 



-1/3 



EL 



kA 1 / 3 . 



(10) 



Notice that the energy in these ridge "wings" or "echoes" 
is negligible compared to the ridge energy which scales as 
kA~ 1//3 This energy is spread over an increasingly 

large area. 

The scaling of the ridge "wings" can be also obtained 
from an extension of the asymptotic analysis of the von 
Karman equations. Ref. |13| determined the scaling of 
the boundary layer solution by rescaling all variables by 
a power of A as in Eq. (|l^) and then requiring that the 
highest derivative terms in Eqs. (|l|) be of the same order 
in A. The exponents of the A-factors that rescaled / and 
y were identical due to the imposed boundary condition 
/ = a\y\ (this a particular way to satisfy Eq. (^|)). Here 
we are seeking the scaling of the large distance solution 
that is pieced together with the ridge solution. Hence, 
we can find such a solution separately for either side of 
the ridge. This allows for arbitrary rescaling factors for 
/ and y since the large distance behavior of / in the case 
of the "wing" is no longer required to be linear in y. The 
matching condition on say d 2 f jdx 2 for some fixed y and 
A — > requires that the "wing" solution for / scale with 
A in the same way as the ridge solution i.e. f = A 1 / 3 /, 
where / is finite in the A — > limit. The same reasoning 
applies to the asymptotic scaling of x- Thus the rescaling 
transformation that is needed to determine the scaling of 
y in the "wing" solution is 



/ = A 1 / 3 /, X = A- 2/3 x, y = \ fJ y, x = A° 3 



(11) 



The longitudinal direction x is not affected by the rescal- 
ing transformation, and the exponent (3 must be negative 
since a positive f3 would reproduce the boundary layer 
scaling. Plugging the scaling ansatz Eq. (jll]) into the 
von Karman Eqs. (Q) we obtain 
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A 1 / 3 



9a; 4 di 2 dy 2 dy A 

= A - 2 /3+l/3-2 /3[ - !/ - ] 



(12a) 



A 2-2/3 



d 4 x 

<9i 4 



2A 



-2/9. 



a 4 x 



dx 2 dy 



-48 ^X 

dy 



TA 



= -\\ 2/ *- 2f) [fJ]. 



(12b) 



Balancing the dominant terms we obtain /3 = —1/3 in 
contrast to (3 = 1/3 for the boundary layer solution jl3) . 
Therefore, the decay length L of the wing strains and 
curvatures scales as VA -1 / 3 in agreement with the pre- 
diction of the energy scaling argument above. The lead- 
ing order behavior of the elastic energy in the wings can 
be found by substituting the rescaled variables into Eqs. 
(0). We obtain i5 W i ng ~ kA 1 / 3 in accord with the energy 
scaling argument. In Sec. VI we numerically verify the 
existence of the long range decay of the curvatures and 
stresses away from the ridge. The decay length is found 
to scale as predicted in the small thickness limit. 



V. RIDGE UNDER EXTERNAL FORCING 

To assess the relevancy of the results obtained for the 
"minimal ridge" to ridges in a crumpled membrane, we 
must first discuss ways in which the effective bound- 
ary conditions for a ridge in a crumpled sheet differ 
form that of the "minimal ridge," and second, we ought 
to determine how these differences affect such relevant 
ridge properties as the coefficients in the thickness scaling 
laws. The additions and changes to the rectilinear frame 
boundary conditions for the "minimal ridge," that distin- 
guish them from realistic boundary conditions for a ridge 
in a crumpled sheet, can consist of a) stresses applied at 
the boundary, b) torques applied at the boundary, and c) 
distributed normal forces that arise when distant parts 
of the crumpled sheet press on the ridge. 

In this section we discuss the linear response of the 
"minimal ridge" to external perturbations. A priori, 
there is no reason to believe that the asymptotic scal- 
ing of the linear response moduli is independent of the 
details of the boundary conditions that create and, more 
importantly, maintain the ridge in the process of exter- 
nal loading. In other words, if we have chosen to main- 
tain constant normal forces on the boundaries instead of 
maintaining constant curvature, the scaling of the linear 
response moduli with the thickness could be different. In 
fact, numerical evidence suggests that the scaling of the 
ridge stiffness with respect to compression does indeed 
depend on the way boundary conditions are maintained 
as the ridge is being distorted. One must therefore ad- 
dress the applicability of the results derived in this sec- 
tion to determination of the elastic response of crumpled 
sheets where the details of the boundary conditions are 



not known. We are unable to address this question here 
beyond suggesting that the linear response of a regular 
tetrahedron considered in Sec. VII might reflect the situa- 
tion in the crumpled sheet given our view of the crumpled 
membrane as a collection of vertices, ridges and facets. 



A. Linear response to specified boundary forces 

We present a treatment that allows one to construct a 
consistent perturbation expansion in the small external 
forces that act at the boundary. The case of the dis- 
tributed external forces can treated in a similar manner 
that will be discussed below. This method of treating 
perturbations is by no means unique. It allows one, how- 
ever, to easily establish asymptotic scaling of the energy 
correction in the limit of the vanishing membrane thick- 
ness. The essential idea is to reformulate the problem in 
terms of some new functions fi and \i that are subject 
the unchanged "minimal ridge" boundary conditions Bq 
but satisfy modified equations. Consider the case when 
the full boundary forces B ext (including the applied ex- 
ternal forces) can be decomposed into B cxt — B + SB. 
Let f e and \e be the solution to the boundary value prob- 
lem SB that includes only the small additional forces (i.e. 
the strip is not bent). We then seek the solution to Eqs. 
(jj]) subject to the full boundary conditions in the form 



/ = fi + fe 



X = Xi+ Xe 



(13) 



so that fi and Xi are subject to the same boundary con- 
ditions Bq as the undisturbed ridge solution but satisfy 
modified equations 



V 4 fi=[Xi,fi} + \Xe,fi} + \Xhfe] 



A 2 V 4 Xl 



" 2 [/«! fi\ [fe: fi_ 



(14a) 



(14b) 



For a sufficiently small perturbation one can construct a 
series expansion around the unperturbed solution. 

The equations for the first order corrections to the 
ridge solution are linear and inhomogeneous. The coeffi- 
cients as well as the source terms are proportional to the 
second derivatives of the unperturbed "minimal ridge" 
solution. These second derivatives are proportional to 
the stresses and curvatures in the unperturbed ridge. 
According to of Ref. |l3| the dominant (dimensionless) 
stresses and curvatures are of order A -1 / 3 in the region of 
the ridge of width A 1 / 3 . Sec. IV of this article determines 
the behavior of the stresses and curvatures in the rest of 
the sheet. The dominant non-dimensional stresses and 
curvatures in the ridge "wings" are of order A 1 / 3 . These 
"echo" disturbances are spread over a large region of size 
L ~ A -1 / 3 so that the elastic energy in the "wings" is 
negligible compared to the ridge energy. The correction 
to the ridge solution will possess the same qualitative 
features as the source terms in the equations (|J) that 
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determine them. Since the behavior of the source terms 
in Eqs. (|lj) is completely determined by the "minimal 
ridge" solution we postulate that the perturbing energy 
is confined to the region of the ridge of width A 1//3 . We 
will use this feature later to determine the appropriate 
integration domain for the ridge energy correction due to 
the external forces. 

At this point we must mention how distributed forces 
could be included in this treatment. On the one hand, we 
could find the effect of the external forces on the flat strip 
with free boundaries and then proceed with the deriva- 
tion as above. On the other hand, distributed normal 
forces can be directly added to the first von Karman 
equation since it states the normal force balance on an in- 
finitesimal element of the sheet dxdy. External in-plane 
forces lead to a redefinition of the Airy stress function 
X- It may not be at all possible to define a stress func- 
tion when in-plane external forces are present, however. 
External distributed torques result in a redefinition of / 
if its definition is indeed possible. This is certainly the 
case for small deflections when Monge coordinates can 
be used @. 

The method presented above of splitting the bound- 
ary value problem into two is applicable when several 
conditions are satisfied. First, note that if the equations 
were linear, the result of such splitting would be trivial 
since the equations for fi and Xi would be the same as 
the original equations and the result of the splitting the 
problem into two amounts to the principle of superpo- 
sition. Second, the perturbed boundary conditions B ext 
that include the external force must be compatible with 
the boundary conditions Bq. For compatibility we re- 
quire that both sets of boundary conditions are supplied 
in terms of the same functions (and derivatives) of / and 
X at the boundary. Thus the change from Bq to B ext 
amounts to a change in the boundary values of the spec- 
ified functions. We denote the change in these boundary 
values schematically as SB = B ext — Bq. This condi- 
tion assures that the boundary value problem SB for f e 
and Xe is well-posed. Third, we assumed that the part 
of the boundary conditions Bq that create the ridge is 
unchanged in the process of loading. 

It is unknown to us whether these conditions are sat- 
isfied in a crumpled sheet. We will see below that the 
prediction based on this method, for the linear response 
to point forces acting at the vertices of a regular tetrahe- 
dron, is incorrect. A possible explanation is that precisely 
the condition for Bq to be held constant in the process 
of loading is violated in the tetrahedron. The effective 
boundary conditions for each ridge in a tetrahedron, the 
stresses in the facets, for example, may change in propor- 
tion to the forces acting at the vertices. However, it might 
by possible to modify the conclusions of this method for 
a general case of changing boundary conditions. Certain 
features that emerge from the perturbation scheme are 
likely to survive for a general loading. In particular, the 
effect of the perturbation is likely to be confined to the 
ridge region as above. 



Another limitation of this method stems from the ar- 
bitrariness of the choice of the external load. In other 
words, of all possible ways to load the ridge, the one in 
the ridge's "weakest direction" is relevant to determina- 
tion of the elastic response of crumpled sheets. It may 
happen that due to high symmetry, the chosen load does 
not "have a component in the weakest direction" of the 
ridge. As a result, the stiffness of the ridge in response to 
such a load is qualitatively greater than its stiffness with 
respect to a generic load. The concepts in quotes above 
can be made rigorous. We discuss various types of load- 
ing in a later subsection with regard to their relevance to 
the determination of the weakest ridge modulus. 

B. Point forces applied at the ridge vertices 




FIG. 2. A long strip of width X is bent by normal 
boundary forces (one quarter of which are shown) and is 
compressed by forces F applied at the vertices. 

Analysis can be carried further when the external in- 
fluence is characterized by a single small dimensionless 
parameter F. In particular, for the most part of the fol- 
lowing development, we explicitly consider the case of 
compression of the ridge by a pair of point forces acting 
at the ridge vertices as shown on Fig. 2. Presumably, 
there are loads for which the asymptotic A-scaling of the 
solution to the "flat" problem f e and Xe is different from 
that in the case of the point forces. This fact however 
introduces but a few changes in the following derivation. 

A flat strip compressed by a pair of point forces acting 
at the boundary remains flat for small enough forces. In 
addition, the stress function is linear in the applied force 
F 

f e (x,y)=0 Xe(x, y) — FX~ 1 4>(x, y), (15) 

where the function (j)(x, y) depends neither on the force 
F nor on the dimensionless thickness A. Above a buck- 
ling threshold for the compressing forces, the flat solution 
Eq. ( |i"5| ) becomes unstable to small perturbations. An- 
other buckled solution exists and may have a different 
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dependence on the force F. This fact does not affect the 
linear response analysis, however, since the flat solution 
is unique for small F. 

We proceed with the boundary layer analysis following 
Rcf . E| . This involves rescaling all variables by a power 
of A 

fi = X fi, Xi = X 5 Xi , x = X°x, y = X fJ y 7 (16) 

with (5 = 2/3 and f3 = —1/3 chosen in such a way as to 
balance the powers of A in front of the terms in the Eqs. 
(|l4f ) that are relevant in the boundary layer. As in Ref. 
|p3| , we seek the solution to the rescaled equations as a 
series expansion in A 2 / 3 



/ i = /f ) +A 2 / 3 /i 1) + A 4 / 3 /f ) + 
h = x f + X 2 ^ + A 4 /3 x f ) + . 



(17a) 
(17b) 



/( o) = ; (o) + fe „ A i/ 3(/o + p h +P 2 h) (2Qa) 

x (o) = x (0) + Xp ~ A~ 2 / 3 (xo + F X i + F 2 % 2 ) + A- 1 ^. 

(20b) 

The terms linear in F completely characterize the lin- 
ear response of the ridge to this particular type of load- 
ing. In general, according to the postulated confinement 
property of the perturbing energy, second derivatives of 
fi and xi possess the same qualitative features as those 
of /o and Xo- I n particular, the dominant second deriva- 
tives of fi and xi are significant in the ridge region of 
characteristic width A 1 / 3 . 

Substituting the F-expansions into Eqs. (||) we obtain 
the expressions for the bending and stretching energies 
accurate up to terms of order F 2 



The equations for the zeroth order terms and Xi°^ 



E° + E°F 



b 



read 



-Estr — Eq + E{F - 



E° 2 F 



E^F 



(21a) 



(21b) 
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j-^(o) y(o)j _j_ xi/Sp^j-i ) (18a) ^he ener gy °f t ne unperturbed ridge as in Ref. |T 



El 



dx 



(0) 



7 A 



(18b) 



Eq ~ A 1 / 3 . Other E^s are sums of integrals that 
involve second derivatives of fi, Xi and 0- For example 



Let us examine the source term in detail. 



(o) 
i >• 



d 2 jf > d^(i,y) A _ 2/3 a 2 /i 0) d*cf>(x,y) 



F / dxdy 



A -l/3+2p 



d 2 Xi 
dy 2 



d 2 Xi 
dy 2 



dx 1 dy 2 



dy 



it" 



dx 2 



2A- 1 / 3 . 



dxdy dxdy 



(19) 



The derivatives of <f> are evaluated at x = x and y — 
A 1 / 3 ?/. Since the rescaled variables are finite in the A — > 
limit, 0's behavior near y — determines the leading or- 
der behavior of its derivatives in Eq. ( |l9|) in the A — > 
limit. The magnitude of the perturbation term is thus 
governed by F = X P F where the value of the exponent 
p is determined by the behavior of the derivatives of <p 
near the ridge. For example, suppose all derivatives of <j> 
are finite at y = 0. This is indeed the case for the com- 
pression of a strip by a pair of point forces Jl8[] . Then, 
the second term in Eq. (|l^) dominates so that p = —1/3. 
This behavior need not be generic, since the derivatives 
of <p may vanish at y = by reasons of symmetry. Then 
the second term in Eq. (|l9|) may vanish so that a higher 
order term in A dominates. This results in more positive 
value of p. 

The solution f- ' and x| ' to Eqs. ( |ilj| ) can be sought 
as a series expansion in the small parameter F. Thus, 
the zeroth order term in the A 2 / 3 -expansion of the so- 
lution to the von Karman equations in the presence of 
the external forces acting on the ridge is given by (up to 
quadratic terms in F) 



+ F 2 j dxdy {V 2 4>f. (22) 

We must remark that the terms quadratic in F that in- 
volve f 2 and xi cancel each other so that the linear in 
F correction to the potentials / and x determines the 
quadratic correction to the total energy. This situation 
is common to all linear response problems. 

The leading order behavior of E± and E% in the small 
thickness limit can be found by assuming, as before, that 
the corrections to the ridge solution /i and xi are con- 
fined to the ridge region. The domain of integration in y 
in the integrals involving these ridge corrections as in Eq. 
(^||) is thus of order A 1 / 3 . We can therefore obtain the 
asymptotic scaling of the coefficients in Eqs. (^T|). The 
bending and the stretching pieces scale with A in the same 
way. The correction to the total energy is E 2 F 2 



E (b,s) ^ A (M) A -l/3+p + j g( 



(b,s) 



E<> — Eo 



A 2 \- 1 / 3 + 2 p + B 2 X p + C 2 



(23a) 
(23b) 



where A\ etc. are arbitrary coefficients that do not de- 
pend on A in the A — > limit. If external stresses and 
curvatures given by the derivatives of Xe and f e respec- 
tively are zero on the ridge line and do not increase suffi- 
ciently rapidly away from the ridge, the exponent p will 
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be positive and large so that the energy correction will be 
dominated by the last terms in the expressions for E\ and 
E 2 . These terms are identical to the energy of a flat sheet 
acted on by the external forces. If p < 1/6, however, the 
energy change is dominated in the small thickness limit 
by the interaction of the ridge with the external forces 
given by the corrections to the ridge solution /i, xi- 

Eqs. (^3|) predict that both kinds of energy separately 
depend linearly on F. We demonstrate numerically in 
Sec. VII that upon application of point forces at the ridge 
vertices, each energy correction does indeed depend lin- 
early on the applied force with the coefficient that scales 
with A as predicted. The total energy can depend only 
on the square of the applied force since the boundary 
stresses and torques vanish for the undisturbed ridge so 
that the work done by the external force is exactly the 
change in the ridge's elastic energy. This means that the 
linear terms in F in the expressions ( p3| ) for -Ebond and 
i? s tr must cancel each other, i.e. Ef — —E\ to all orders 
in the A 2 / 3 -expansion {A\ = -A\, Bf = -B\ etc.). The 
numerics reported in Sec. VII convincingly show that it 
is indeed the case. The A-scaling of the ridge's elastic 
constant E2 given in the Eq. (|23|b) has also been shown 
to be consistent with the numerics. 

When external forces distort the ridge, additional lon- 
gitudinal strain results. We first define a "vertex-to- 
vertex" strain j v as the amount by which the vertices 
move closer together divided by X. We compare this 
quantity with the additional strain in the sheet that is 
the change in j xx ~ (l/Y)d 2 x/dy 2 . It can be found from 
the first order in F correction to \ Eq. (p(i|). We can in- 
fer 7„ by inspection of the energy- force relation, since the 
work done by F is the change of the energy E 2 F 2 . 

lv ~ \E 2 F ~ {A 2 X 2 / 3+2 p + B 2 X 1+ p + C 2 X)F = G^F, 

(24) 



whereas from Eq. fl20| ) we obtain the change in the lon- 
gitudinal strain in the sheet 



add 

1 xx 



X 2/3+pp 



(25) 



Eq. ( pq ) gives another way to estimate the exponent p 
numerically. 

This way of looking at the problem motivates the 
following consistency check of the perturbing energy 
confinement. An additional way of determining the 
A-scaling of the total energy correction via E 2 F 2 ~ 
G(lxx d ) 2w X/n ~ A~ 2 / 3+2p w must reproduce the scal- 
ing of E 2 obtained by explicit integration as in Eq. (23). 
Here w = w/X is the characteristic width of the region 
in which the perturbing energy is confined. Calculated 
by this method, the width id ~ A 1 / 3 is independent of 
the exponent p and scales with A in the same way as the 
ridge width. 



C. Buckling threshold 

An important property of the ridge solution that is in- 
accessible by this simple-minded perturbation scheme is 
the force required to buckle the ridge. In principle, one 
must solve Eqs. (|l8|) and then perform a linear stability 
analysis of the solution to determine when the ridge will 
buckle. This task is intractable analytically due to the 
complexity of the equations. The asymptotic scaling of 
the buckling threshold may nevertheless be anticipated 
using the following argument. The undisturbed ridge so- 
lution /o and xo is linearly stable against shape pertur- 
bations. Therefore, changes in this solution that destroy 
stability are likely to be of the order of the solution itself. 
To induce such changes, the additional terms in the Eqs. 
( |l8| ) for the loaded ridge must be comparable in magni- 
tude to the rest of the terms in the equations. Hence, 
the parameter F that controls the magnitude of the ad- 
ditional terms in the Eqs. (|l^) has to be of order unity 
F ru 1. This means that in the case of ridge compressed 
by vertex forces, the buckling threshold force scales as 



F n , 



A" 



(26) 



This argument has a serious flaw. The scaling of the 
"flat-problem" solution f e and Xe may be different for 
large forces F since additional solutions to the "flat- 
problem" may appear. This may result in a different 
value of the exponent p and thus a different scaling of the 
buckling threshold. The conjectured scaling of the criti- 
cal load has two important implications. First, the addi- 
tional longitudinal sheet strain at the buckling threshold 
given by the Eq. (|2^) is of order of the strain that ex- 
isted in the ridge prior to compression. Second, the ridge 
energy correction at the buckling threshold value of the 
external force has the same asymptotic thickness scaling 
as the energy in the undisturbed ridge. This conclusion 
is supported numerically in Sec. VII. A corollary of this 
statement is that the applied load on a ridge does not 
have a decisive effect on its energy. The ridge buckles 
when its energy changes by a finite fraction. We will 
see in the following subsection that this argument ap- 
plies to other types of loading as well. Therefore, any 
two unbuckled ridges in a crumpled sheet with compara- 
ble lengths and dihedral angles should have comparable 
energies, even though they have different loads. 



D. Other types of loading 

To describe elastic properties of the ridge in a way that 
is relevant to determination of the structure of the crum- 
pled sheet we must investigate other types of loading. 
We first consider other loads that can be treated with 
the perturbation scheme developed above. These include 
normal forces acting on the "minimal ridge." Having 
done that, we anticipate that in a crumpled sheet the 
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assumptions that lead to the scaling laws for the energy 
corrections are violated. 

Let us first use the perturbation scheme developed 
above to discuss other loads. The goal is to determine 
the "weakest modulus" of the ridge alluded to above. 
The perturbation treatment of this section relates the 
linear response of the flat sheet to the linear response 
of the ridge. We therefore anticipate that any pertur- 
bation that causes a large response in a flat sheet will 
also be relevant to the determination of the "weakest" 
ridge modulus. Since the bending rigidity k ~ Yh 3 of 
an elastic sheet vanishes faster in the h — > limit than 
the stretching modulus G ~ Yh, any perturbation that 
causes the sheet to bend will create a large response. We 
therefore consider normal forces on a ridge presented on 
Fig. 3. 




FIG. 3. Same strip is bent by normal forces applied at 
the vertices as well as the middle of the ridge. 

Let us first describe the effect of these forces on a flat 
sheet. The scaling of the curvature potential f e with 
force F can be found by noting that the torque of the 
external force is balanced by the flexural moment of the 
sheet due to its curvature kX d 2 f e /dx 2 ~ FX. There- 
fore, f e = \~ 1 Fip is controlled by a single dimensionless 
parameter \~ X F. Since the sheet is bent in the short di- 
rection, we can use the arguments of Sec. IV to show 
that this curvature persists up to a distance of order 
L ~ XX^ 1 / 3 . The treatment of Sec. IV also determines 
the stresses in the strip. The dominant component of 
the stress is a yy ~ Yh(a/L) 2 where a ~ f e ~ X\~ 1 F. 
Hence, the stress is of the second order in the applied 
force F. The modified von Karman equations ( |l4|) will 
thus have a first order in F source term due to f e and a 
second order in F source term due to % e . One can now 
repeat the steps that lead from the modified von Karman 
equations ( |l4| ) via a rescaling transformation Eq. (|l^) to 
the perturbation expansion Eqs. (21). The only differ- 



ence is the form of the flat solution f e and \e- 

Together with knowledge of the small y behavior of 



the function ip one can determine the exponent p which 
controls the ridge stiffness under this particular type of 
loading. The important second derivative d 2 (p/dx 2 is fi- 
nite at the ridge. Therefore, p = —4/3. This means first 
that the linear response of the ridge to normal forces is 
qualitatively greater than to in-plane forces. Second, the 
normal force required to buckle the ridge F cr u ~ A 4 / 3 
is much smaller than in the case of longitudinal com- 
pression. The conjecture that the ridge buckles when its 
energy changes by a finite fraction can be made on the 
same basis as for the case of the longitudinal compression 
of the ridge. 

We have thus shown that depending on the symme- 
try of the applied load, the linear response of the ridge 
can vary considerably. It is unknown to us at this point 
whether the resistance of a crumpled sheet to further 
compression is determined solely by the "weakest" ridge 
modulus. It is certainly not inconceivable to imagine a 
situation in which a number of ridges in a crumpled sheet 
that form a kind of a "structural skeleton" are loaded 
in such a way that their linear response is given by a 
stronger modulus. 

Let us finally discuss the applicability of the pertur- 
bation scheme developed here to ridges in a crumpled 
sheet. We anticipate the boundary stresses for the ridge 
in a crumpled sheet to change by an amount proportional 
to a load. We cannot, therefore, utilize the perturbation 
method, developed in this section, to determine the scal- 
ing of the energy correction and the buckling threshold. 
Certain features of the conclusions that emerged from 
the perturbation analysis are likely to persist, however. 
We anticipate that, the perturbing energy is confined to 
the ridge region. In addition, the change in the ridge 
longitudinal strain (the dominant ridge strain) must be 
proportional to the movement of the sheet caused by the 
external forces. The "vertex-to- vertex" strain is an exam- 
ple of such movement. The coefficient of proportionality 
is likely to scale as a power of the thickness A. In most 
situations, one can make a reasonable guess to what that 
power is. For example, the numerics reported in Sec. VII 
show that when a tetrahedron is compressed by point 
forces applied to its vertices, the additional ridge strain 
scales in the same way with A as the "vertex-to-vertex" 
strain. This assumption, together with the notion of the 
perturbing energy confinement, is sufficient to determine 
the scaling of the energy correction. 



VI. RIDGE INTERACTION 

Ridges in a crumpled elastic sheet can interact in two 
distinct ways. First, distant parts of the sheet can press 
against a ridge. This type of interaction can be discussed 
using the framework developed in Sec. V. However, a 
quantitative understanding of the distribution of forces 
through the crumpled sheet is needed. Second, nearest 
neighbor ridges can influence one another through the 
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strains and curvatures they create in the sheet. We will 
concentrate on this particular aspect of the ridge interac- 
tion since we are able to obtain quantitative conclusions. 

Two length scales characterize a ridge. Most of the 
ridge's elastic energy is confined within a strip of width 
w ~ XX 1 / 3 around the ridge. Small curvatures and 
strains persist up to a distance L ~ XX^ 1 ^ 3 from the 
ridge. In a crumpled elastic sheet typical distances be- 
tween ridges D are of the order of their length: D ~ X . 
Assuming that all ridges have the same characteristic size 
X we are led to a conclusion that w <C D <C L for most 
ridges. Therefore, the small residual strains and curva- 
tures of a ridge will influence its neighbors in a way that 
is calculable within the framework of Sec. V. In fact, in 
Sec. IV we argued that there is a small transverse strain 
lyy ~ fl2 1^ 1 ~ ^ 3 present in the sheet. This strain can 
be thought of as resulting from an external force F ~ A 1 / 3 
that is stretching the sheet in the direction transverse to 
the ridge. Therefore, the stress potential function <fi has 
a finite derivative d 2 <fi/dx 2 at the ridge. According to 
the prescription of Sec. V, p — —1/3, so that the energy 
correction scales as 



SE ~ kX-^+^F 2 ~ kX- 1 ' 3 . 



(27) 



Note that the energy correction scales the same way with 
A as the ridge energy. This means that even though the 
ridge "wings" carry a negligible amount of energy, they 
cause the total energy of the system to be changed by 
a finite fraction. This is not surprising since the trans- 
verse "echo" stresses discussed in Sec. IV are of the same 
magnitude as the transverse stresses in the ridge. The 
interaction energy of ridges in a strip geometry is stud- 
ied numerically in the following section. It is found to be 
a few percent of the ridge energy in the small thickness 
limit. 



VII. NUMERICAL SIMULATION OF THE 
RIDGE PROPERTIES 

In this section we use a lattice model of an elastic sheet 
of Ref. 0. This model was used by the authors and 
others to verify a number of the ridge scaling properties 
|p^ , ^3| . The sheet is modeled by a triangular lattice of 
springs of equilibrium length 6 and spring constant K . A 
bending energy J(l — ni ■ n 2 ) is assigned to each pair of 
adjacent lattice triangles with normals rii and fi2. When 
the strains are small compared to unity and radii of cur- 
vature are large compared to the lattice spacing 6, this 
model bends and st retches like an elastic sheet of thick- 
ness h — b^SJ/K and bending modulus k = JV3/2. 
Ridges could be created either by imposing appropriate 
boundary conditions to a long strip of the simulated ma- 
terial or by introducing disclinations. We studied both 
types of shapes. First, we applied forces to the particles 
located on the long boundaries of a strip so as to con- 
strain them to lie in different planes on each side of the 



ridge. The angle between the normals to these planes is 
2a. Second, we connected a triangular piece of this sim- 
ulated material into a regular tetrahedron so that each 
edge then became a ridge. A sequence of minimum en- 
ergy shapes of different dimensionless thickness A were 
obtained with the use of a conjugate gradient routine. 
External forces were then applied to the ridge vertices 
and linear response to compression was measured as well 
as the buckling characteristics. 




FIG. 4. Ratio of the total bending to the total stretch- 
ing energy vs. the dimensionless thickness A for a regular 
tetrahedron of edge length 1006. The asymptotic limit of 
5 is approached for A ~ 10 ~ 4 . 



We first tested the virial theorem that predicts that if 
all of the elastic energy is concentrated in the ridges then 
the total bending energy is five times the total bending 
energy. In Fig. 4 we plot the ratio of the total bend- 
ing to total stretching energies for the tetrahedron with 
edge length of 1006. This ratio approaches the predicted 
value of £'bcnd/-E'str = 5 within a few percent for the val- 
ues of dimensionless thickness A < 10~ 4 . For the smallest 
values of A on the plot the ratio deviates from the asymp- 
totic value due to lattice effects. Thus the virial theorem 
is obeyed within our measurement precision. 

The predictions of the scaling behavior of the ridge un- 
der external loading was done using a 43.36 by 1006 strip 
bent by boundary forces as well as a tetrahedron with an 
edge length of 506. We applied point forces to the ver- 
tex particles in either case in such a way as to compress 
the ridge. We first compressed a flat sheet to verify that 
stresses are finite on the line y = 0. Then, the prediction 
for the exponent p defined in Sec. V is p = —1/3. The 
correction to the total elastic energy was found to vary 
on F 2 as anticipated. The coefficient E<z of F 2 should 
scale with A according to Eq. ( |23| ) . Fig. 5 is a plot of E 2 
as extracted by a quadratic fit of the correction of the 
total elastic energy per ridge for the tetrahedron (trian- 
gles) and the strip (squares) as a function of A. The inset 
of Fig. 5 shows the fit to the quadratic dependence of the 
total energy on the applied force. 
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FIG. 5. Thickness dependence of the coefficient E% in 
E/n ~ E0 + E2F 2 for the strip of dimensions 436 by 1006 
(squares) and tetrahedron of edge length 506 (triangles) . 
The solid lines are fits according to the Eq. J23]) with 
p = —1/3 for the strip and p = for the tetrahedron. 
Inset: Plot of the total energy of a strip as a function of 
the applied force F. Solid line is the quadratic fit to the 
first five points. Dashed vertical line marks the buckling 
threshold. 



The numerically accessible range of A does not allow 
for a direct determination of the exponent p from the 
data. However, the strip data are consistent with the 
prediction for p = —1/3. However, the tetrahedron data 
are consistent with p = and inconsistent with p = — 1/3 
when fitted with the scaling form Eqs. (23) and (|24|). A 
given load appears to store qualitatively less energy for 
the tetrahedron than for the strip, i.e. the tetrahedron is 
qualitatively stiffer. This difference might be explained 
by noticing that while in the case of the strip, the bound- 
ary shape Bq that maintained the ridge was fixed under 
the loading, whereas the effective boundary shape for the 
ridges in a tetrahedron changes when forces are applied. 
In addition, tangential stresses act on the the effective 
ridge boundary in a tetrahedron, whereas only normal 
boundary forces are present in the strip geometry. There- 
fore, the decomposition method of Sec. V cannot be used 
for the tetrahedron. The scaling of the tetrahedron stiff- 
ness implies that the additional strain due to the action of 
the compressive force is confined to the ridge region and 
scales the same way with the thickness A as the "vertex- 
to-vertex" strain. These results are different from the sit- 
uation in the strip. There, the "vertex-to- vertex" strain 
may be relaxed through bending so that the additional 
strain in the sheet is qualitatively weaker than the move- 
ment of the vertices would dictate if such relaxation was 
not possible. 



FIG. 6. The coefficient E\ of the linear dependence 
of the total bending energy on the force F for the 436 
by 1006 strip (squares) and the tetrahedron of size 506 
(triangles) . The solid lines are again fits according to the 
Eq. ( |23| ) with the same values of the exponent p as in 
Fig. 5. 

In Fig. 6 we plot E\ for both shapes extracted by a 
linear fit of the total bending energy as a function of F. 
The data are consistent with the prediction Eq. ( p3|) for 
the same values of the exponents p as obtained from the 
scaling of E^- Finally, we graph the coefficient G -1 of 
the linear dependence of the induced "vertex-to-vertex" 
strain on the applied force in Fig. 7. The data are con- 
sistent with the prediction Eq. (|2J). For comparison let 
us consider a flat strip. Its dimensionless inverse stiffness 
G _1 to compression diverges as A -1 , whereas its stiffness 
to bending vanishes as A. 



0.0025 




0.0004 



FIG. 7. The coefficient G^ go of the linear F depen- 
dence of the "vertex-to- vertex" strain created by the ex- 
ternal compressive force F for the strip (squares) and the 
tetrahedron (triangles). Solid lines are fits according the 
Eq. (p4|). The tetrahedron ridge stiffness G r id ge is quali- 
tatively greater in the limit of the vanishing thickness. 

An important question that needs to be settled numer- 
ically is the buckling threshold of the ridge. The vertex 
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forces were increased until the shape underwent a radical 
change. This seems to imply that buckling is a first or- 
der phenomenon. In other words, there are several stable 
shapes other than the ridge for a range of compressive 
forces below the buckling threshold. Fig. 8 shows two 
such buckled bent strip shapes. A buckled tetrahedron 
shape shown in Fig. 9 exhibits a more complex buckling 
pattern. Shading is proportional to the stretching energy 
density. 




FIG. 8. A roof-shaped strip (left) is collapsed by cor- 
ner forces to obtain two stable buckled shapes as viewed 
obliquely from above. The value of the vertex forces 
must be reduced below the buckling threshold in order to 
achieve stability. Shading is proportional to the stretch- 
ing energy density. 




FIG. 9. A tetrahedron lattice (left) is collapsed by ap- 
plication of vertex forces (right). A complicated buckling 
pattern results. 



In Fig. 10 we plot the ratio of the energy correction to 
the total elastic energy of the undisturbed ridge at the 
bifurcation point of the loaded tetrahedron. This ratio 
seems to be approaching a constant in the A — > limit 
which would agree with the prediction that the energy of 
a ridge can only change by a finite fraction before it buck- 
les. More numerical tests are needed however to establish 
this assertion firmly. We did not attempt to characterize 
the buckling of the tetrahedron in further detail. 
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FIG. 10. Ratio of the energy correction at the buckling 
threshold of the tetrahedron of edge length 506 to the en- 
ergy of the undisturbed tetrahedron as a function of the 
dimensionless thickness A. The empirically drawn solid 
line suggests that the energy correction varies roughly 
linearly with A. 

The large distance behavior of the ridge "wings" was 
tested by using a long 526 by 5006 strip. Fig. 11 dis- 
plays the transverse curvature C yy and the longitudinal 
curvature C xx in the units of X^ 1 along a perpendicular 
bisector of the ridge. A long almost linear decay of the 
C xx is evident whereas C yy decays rapidly to zero. 
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FIG. 11. Transverse ridge curvature C yy (squares) and 
longitudinal curvature C xx (circles) in the units of A -1 
along the perpendicular bisector of the ridge vs. the dis- 
tance from the ridge for a 526 by 5006 strip. 

We extracted the decay length L for a sequence of 
ridges of varying A and plotted the results in Fig. 12 
versus the predicted scaling variable A" 1 / 3 . The linear 
fit shows that the data agree well with the prediction. 
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FIG. 12. Decay length of the longitudinal curvature L 
in the units of X vs. the predicted scaling variable A -1 / 3 
for a 526 by 5006 strip. 

We next investigated the question of ridge interaction. 
Two parallel ridges were created in the strip geometry. 
We varied the distance between the ridges D as well as 
the thickness A. Detailed interaction features that are 
shown in Fig. 13 for the 1006 by 1306 strips for a fixed 
A = 0.0002 depend on the relative orientation of the 
ridges. If they are both concave (or both convex) as 
seen by looking down on the sheet from above (config- 
uration U), they repel at short distances and attract at 
long distances. 

50 I 1 1 1 1 1 




Inter - Ridge Distance D/X 

FIG. 13. Total elastic energy of a two ridge configu- 
ration in the units of k found from a 1006 by 1306 strip 
bent in two places by 90° vs. the distance between the 
ridges D. The thickness is A = 0.0002. Square symbols 
correspond to the ridges that have the same orientation 
(configuration U), whereas the triangles correspond to 
two ridges of the opposite orientation (configuration Z). 

The situation is reversed if the ridges have different 
orientation (one is folded up and the other one down, or 
vice versa: configuration Z). 

The signs of the interactions can be readily understood. 



At short distances, the two interacting ridges become a 
single ridge. If the two ridges have opposite signs, there 
is no ridge when they are brought together and thus the 
total energy is zero. If they have the same sign, the com- 
bined ridge has twice the dihedral angle a as the con- 
stituent ridge. Since ridge energy scales as a 7 / 3 jl3), the 
energy at D = should be 2 7 / 3 /2 times that at D = oo. 
The behavior of the interaction energy at large D is 
also understandable. Here the deformation between the 
ridges is minimal. For opposite sign ridges (configuration 
Z), the curvature fields created in the region between the 
ridges have opposite signs. For same-sign ridges shaped 
like a squared letter U, deformations caused by the two 
ridges at the midpoint reinforce each other. In a medium 
with a quadratic energy functional, such reinforcing de- 
formations lead to a reduction in energy. Conversely, two 
opposite-sign ridges in a Z configuration have opposing 
deformations at the midpoint. This yields a repulsion. 



0.08 




0.0004 0.0008 0.0012 
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FIG. 14. Ratio of the difference of the total energy of 
two ridges and twice the energy of once ridge AE divided 
by the energy of one ridge E as a function of A for a fixed 
inter ridge separation D = 1.15X. The data are obtained 
from simulation of a 526 by 5006 strip bent in two places. 
Squares correspond to the configuration U in which the 
two ridges are have the same orientation, and triangles 
correspond to the ridges of the opposite orientation Z. 

For a complete account of the interaction energy, we 
must consider energy stored not only between the ridges, 
but within each ridge due to the presence of the other. 
To investigate such effects we made a long 526 by 5006 
strip bent by 90° in two places. This creates two parallel 
ridges in the same strip whose relative orientation can be 
changed. In Fig. 14 we plot the difference of the elastic 
energy of the two ridges and twice the energy of one ridge 
AE divided by the energy of one ridge £ as a function A 
for a fixed inter-ridge separation of D = 1.15X. Squares 
correspond to the configuration in which the two ridges 
have the same orientation (U), and triangles correspond 
to the ridges of the opposite orientation (Z). The results 
are consistent with the prediction that this ratio must 
approach a finite constant in the A — ► limit. Therefore 
the prediction that ridge interaction changes the system 
by a finite fraction is convincingly confirmed. 
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VIII. DISCUSSION 

In this article we explored properties of the ridge sin- 
gularity in thin elastic plates that may be relevant to a 
quantitative analysis of crumpled elastic membranes. A 
virial theorem that relates the bending and stretching 
contributions to the total elastic energy has been derived 
from an energy scaling argument and verified numeri- 
cally. The virial theorem affords a useful test of elastic 
energy confinement. When most of the elastic energy is 
confined to the ridges the virial theorem predicts the ra- 
tio of the bending and stretching energies in the small 
thickness limit. 

We have developed a perturbation expansion scheme 
that allows one to calculate the effects of external forces 
on the scaling properties of the ridge singularity. We 
found that the problem can be decomposed into first solv- 
ing for the effect of the external forces on the flat sheet 
Xe a- n d f e an d incorporating this solution into the equa- 
tions that describe the ridge singularity in a way that 
is particularly convenient for a perturbation expansion. 
This method allows one to determine the scaling of the 
energy correction. This correction scales with the applied 
force squared and also a power of thickness that depends 
on the details of \e and fe in the ridge region. Two 
different types of scaling were identified. Generally, an 
imposed strain comparable to the strain in a ridge stores 
an energy comparable to the ridge energy. This leads 
to an effective modulus of order E T - lr i ge /X 3 ~ FA 8 / 3 . A 
weaker modulus is possible for isolated ridges with stress- 
free boundaries since the imposed strain can be relaxed 
in ways not accessible to ridges in a crumpled sheet. 

A feature common to all types of ridge loading is that 
the ridge solution becomes unstable when the external 
forces change the energy of the ridge by an amount that 
is comparable to the original undisturbed ridge energy. 
This discovery gives justification to a claim that the en- 
ergy of a crumpled elastic sheet can be found once the 
ridge network is characterized in terms of the ridge sizes 
Xi and their dihedral angles 9j [13[. 

We have found that the large distance part of the ridge 
solution that matches onto the boundary layer solution 
exhibits scaling with the thickness A. We have estab- 
lished by an energy scaling argument as well as an exten- 
sion of the scaling analysis of the von Karman equations 
that small stresses and strains persist up to a distance 
L ~ AA -1 / 3 away from the ridge. Using the frame- 
work developed in this paper we found that the main 
implication of the ridge "wings" is that ridges located 
at distances D <C L away from each other interact in a 
way that changes their total energy by a small but finite 
fraction. 

In the future we propose to build a model for a quanti- 
tative characterization of the ridge network in crumpled 
sheets that incorporates the properties of ridges uncov- 
ered in this article and makes a prediction for the ridge 
size distribution. We hope to characterize the buckling 
behavior more thoroughly and verify our hypothesis that 



the generic scaling of the modulus is that of the tetrahe- 
dron studied here. 
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